Measurement of the noise spectrum using a multiple-pulse sequence 
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A method is proposed for obtaining the spectrum for noise that causes the phase decoherence of 
a qubit directly from experimentally available data. The method is based on a simple relationship 
between the spectrum and the coherence time of the qubit in the presence of a 7r-pulse sequence. 
The relationship is found to hold for every system of a qubit interacting with the classical-noise, 
bosonic, and spin baths. 
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Introduction. — Implementation of quantum informa- 
tion processing (QIP) requires maintaining the quantum 
coherence of a system during its operation. In reality, 
however, systems are not completely isolated from the 
environment. Interactions with the environment are sys- 
tem noise and cause decoherence. One of the challenges 
in QIP implementation is to control the system appro- 
priately in a noisy environment and/or to reduce the 
noise level by weakening the system-environment cou- 
pling. Strategies for these approaches have been devel- 
oped, including quantum error correction [Il-Qi quantum 
estimation [3,@]j ^'^d dynamical decoupling (DD) [B-17|. 

In DD, the application of multiple 7r-pulses cancels 
out the noise and effectively suppresses the system- 
environment coupling 0. This basic idea of DD comes 
from concepts in pulsed nuclear magnetic resonance 
(NMR) Experiments have demonstrated that de- 

coherence can be suppressed by using methods common 
to pulsed NMR, such as alternating-phase-Carr-Purcell 
(APCP) and Carr-Purcell-Meiboom-GiU (CPMG) meth- 
ods [iqI [20I S ubsequent theoretical 0-23] and experi- 
mental [15l4l7| studies on DD in various kinds of spin- 
and charge-related qubit systems other than NMR com- 
pared several DD methods and suggested that optimiza- 
tion of DD requires knowledge of the noise properties. 

For any strategy (not only DD), knowledge of the noise 
properties is helpful because it can be used to improve the 
strategy. It is therefore important to identify the noise 
properties of the environment. For longitudinal decoher- 
ence (energy relaxation) noise, the noise spectrum can 
be obtained from the longitudinal relaxation time (called 
Ti in NMR) [H, HHlil , and a qubit (two-level system) 
has been proposed as a noise spectrum analyzer [23[. In 
contrast, a method for measuring transverse decoherence 
(pure dephasing) noise has not been established although 
several methods (also using a qubit) have recently been 



proposed [10|, |2J, |25| ■ In Ref. [2^, it was pointed out that 



the spectrum by using the Uhrig pulse sequence [8|. In 
Ref. [25j , a method was described for obtaining the spec- 
trum at the Rabi frequency by using a nearly-continuous 
and on-resonant control field. 

In this Letter we propose a method for measuring the 
dephasing noise spectrum in which a simple sequence of 
equidistant 7r-pulses (such as an APCP or CPMG se- 
quence) is used. The spectrum at frequency 7r/2r, where 
2r is the interval between pulses, is evaluated directly 
from experimentally obtained values of the coherence 
times by using a relationship between the noise spectrum 
and the coherence time for a sequence of a sufficiently 
large number of pulses. We show that this relationship 
holds for the classical-noise, spin-boson, and spin-spin 
bath models. 

Model and noise spectrum. — We use a model of a single 
qubit (spin S — 1/2) interacting with the environmental 
degrees of freedom (bath). The Hamiltonian of the total 
system is given by 



(1) 



the relationship between noise spectrum and coherence 
can be used for estimating the spectrum. In Ref. [lo| . 
a method was described for obtaining the moments of 



where ^ and -ffbath are the bath operators. The Pauli 
matrices of the qubit are denoted by ax, Oy, and Uz. 
This is a pure dephasing model because the interaction 
Hamiltonian ffint = {h/2)^az between the qubit and bath 
includes only ct^. In other words, there is no energy relax- 
ation in this model (Ti is infinite in NMR terminology). 
The spin (qubit) is subject to a static magnetic field V, 
and a random magnetic field (noise) ^ in the z-direction; 
the random field ^ is generated by the bath. 

There are several variations of the model depending 
on the nature of the bath. Here we consider three of 
them. The first is a classical- noise model [13] , where ^ is a 
stationary Gaussian stochastic process (classical random 
variable) with zero mean. (In this case, Hhath does not 
appear.) The second one is a spin-boson model in 
which i = J2j + ^'^d i?bath = J2j where 
bj (bj) is the jth mode annihilation (creation) operator of 
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the bosonic bath, and Xj is the coupUng strength between 
the qubit and the jth mode boson. The third one is a 
spin-spin bath model [2^ in which | = J^j + 
and i?bath = Z]j(V2)cjjS-j, where s\_ — + is^y and 
= — is^y. Here, s^, s^, and s-j are respectively the 
X, y, and 2 components of the Pauli matrices of the jth 
spin in the bath, and fij is the coupling strength between 
the qubit (spin of interest) and the jth spin (in the bath) . 

The noise spectrum is defined as the symmetrized 
power spectral density function of the random field: 



(a) 



dt e^-'^(mm+mm) 



I bath 



(2) 



In the classical- noise model, (•)bath is the average over the 
stochastic variable (^(i) = f (i)). In the spin-boson and 
spin-spin bath models, i{t) = e*-^''-"'*/'^^e~*^''=""*/'^, and 
(•)bath = Trbath(/5blthO average in an equilibrium 

state of the bath (pj^^^i^ = e-'^^^-'V^bath)- 

Pulse sequence and generalized coherence time. — The 
density matrix (which describes the state) of the qubit is 
denoted by ps- In the classical- noise model, ps is given 
by the average of ps over the stochastic variable: ps — 
('°'*)bath' ^h^'"^ Ps is the qubit state with one realization 
of the stochastic variable. In the spin-boson and spin- 
spin bath models, ps is given by the partial trace of the 
density matrix p of the total system: ps = Trbath/5- 

The coherence is quantified using a non-diagonal com- 
ponent of the density matrix of the qubit: ps^_| — 

(+|/5s|— ). Here, |±) is the eigenstate of ct^ corresponding 
to the eigenvalue ±1. Experimentally, the coherence is 
measured using the transverse (a; and y) components of 

the qubit (spin) and the relations (ax) = 2Yieps^-\ and 

i^y) = -2Ini/3s,+-- 

Now we consider a situation in which ideal (instan- 
taneous) TT-pulses (about the x or y axis) are repeatedly 
applied to the qubit. As shown in Fig.[IJa), the pulses are 
applied at times ti,t2, ■■■,tm where n is the total number 
of pulses. After the application of the sequential pulses, 
we measure the normalized coherence at time t (> t„): 
Wit) = |p,,+_(t)|/|p.,+_(0)l. 

When we apply a sufficiently large number of pulses 
(keeping the minimal inter-pulse time fixed, so that t is 
also sufficiently large), the coherence exhibits an expo- 
nential decay (as shown in the models later): 



Wit) ^ exp (-t/Ti^) 



(3) 



This time dependence enables us to define uniquely the 
coherence time T2 for a multiple-pulse sequence, which 
we call "generalized" coherence time (for a reason men- 
tioned later). T2 differs from conventional coherence 
time T^^, which is defined using spin echo (SE) exper- 
iments (T^ > in most cases). Generally, it can be 
shown that depends on the pulse sequence. 

In the three models, the time evolution of the normal- 
ized coherence in the presence of a pulse sequence can be 



, 

?! ^2 h h h 



1 



I I 

H 1- 



I I 
H ^—^t 



FIG. 1. (a) Illustration of a 7r-pulse sequence, (b) Function 
ft{t') corresponding to illustration (n: even). 



expressed as (weak coupling conditio n /i, - ^ ujj_is neces- 
sary for the spin-spin bath model) 0- 13 26, 27| 



Wit) = exp £ ^Siu^)\Uu^)\^'^ , (4) 

/OO 
di'e-*7t(i'), (5) 
-00 

n 

Mt') = - *')e(t' - tk), (6) 

k=0 

where tg = 0, tn+i = t, and {tk}^=i is the set of times 
when the pulses are applied, and 6(T) is the Heaviside 
step function. As shown in Fig. [Ijb), ftit') depends on 
the pulse sequence and takes a non-zero value (+1 or — 1) 
only for <t' <t. 

Relationship between Tf^ and S. — Before present- 
ing the main results, we show the relationship be- 
tween Tf^ and S. For the SE pulse sequence (7r/2(a;)- 
T-7r(a;)-T-signal), we can easily find that |/t(a;)p — 
(16/a;^) sin''(a;r/2). Substituting this into Eq. (jH) and 
using the asymptotic behavior (l/w^) sin'*(ti;r/2) ~ 
(7r/4)(5(a;)r as t — > cxd, we get the asymptotic r- 
dependence of the SE coherence VF^^(2r), 



VF^^(2t) - exp 



-S'(0)2r 



00, 



which yields the relationship between and S: 



(7) 



(8) 



Although this formula is common in the field of NMR 
[lil [iij l. it has a significant implication; we should use 
the asymptotic exponential decay of the SE coherence 
to define Tf'^, as seen in Eq. dH) [2?, "SO^ . This is con- 
trastive to the usual evaluation of coherence time in SE 
experiments, in which the behaviors of the SE coherence 
at smaller r are generally used. Their use would not pro- 
vide a unique definition of coherence time because the 
functional form of the SE coherence at smaller r depends 
on the systems. [In some systems exponential decays are 
observed with certain time constants (different from ^ 
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in some cases), while in other systems Gaussian decays 
are observed with certain time constants (different from 
T|E).] 

Relationship between T2 and S. — In the long time 

limit (large n limit) keeping the minimal inter-pulse time 

fixed, we get the power spectral density function of ft{t')- 

I ~ 1 2 I ~ 1 2 

I{uj) EE limt^oo(lA)|/i(w)| . Hence |/t(w)| ^ I{uj)t as 

t 00. From this and Eq. (U), we find that the asymp- 
totic behavior of the coherence is an exponential decay: 



W{t) ~ exp 



^SiL.)I{L0)t 
2,71 



as i — >■ 00. 



(9) 



The necessary number of pulses to extract this asymp- 
totic exponential decay is independent of the noise spec- 
trum because ftit') is determined only by the pulse se- 
quence. Comparing this equation with the definition of 
T2 [Eq. ©], we obtain 



do; 
2^ 



S{u)I{uj). 



(10) 



If the sequence consists of equidistant pulses with 
tk+i - tk ^ 2t {k ^ l,2,...,n - 1), ft{t') is a periodic 
function of t' with a period of 4t (for < t' < t). In 
this case, the power spectral density function I{uj) (with 
fixed r) becomes 



I{uj) = 27r 



E 



(11) 



where cUm = vnir/iT, and Cm is the Fourier coefficient 
oi ft{t'): C,„ = (1/4t) /o^'di' e'^^'' ft{t'). Substituting 
Eq. (dll) into Eq. dTU]), we obtain 



^ 00 



(12) 



m=0 



As an example, we investigate the use of the APCP 
pulse sequence: 7r/2(a;)-{T-7r(a;)-2T-7r(a;)-T}"/^-signal (n: 
even). For this pulse sequence. 



1 for AkT <t' < {4k + 1)t 

and {4k + 3)r < t' < 4{k + l)r, 
-1 for (4fc -I- l)r < < (4fc + 3)t, 
otherwise. 



where k — 0,l,2,...,n — 1. This yields |Cm| = 
{4/TT^m^)6m,2i+i {I — 0, 1, 2, ...). Hence we finally obtain 



1 4 
7^ " 



1 



n^-{2l + iy 



;S{u:2i 



(13) 



Note that we obtain the same result for the Carr-Purcell 
and CPMG sequences because we assume ideal (instanta- 
neous) pulses. [To avoid pulse-error accumulation in ac- 
tual experiments (with non-ideal pulses), one should use 



CPMG or APCP.] Since the factor 1/(2Z + if is smaller 
for larger I, we can approximate the above equation into 



1 

T'L 
^2 



;S{7t/2t), 



(14) 



if S{uj) rapidly decreases as uj increases. Equations (|T3|) 
and (I14|) are the main results of this Letter. 

These relationships are qualitatively explained as fol- 
lows. In many systems, coherence time is dominated by 
the lowest-frequency component of the noise spectrum 
(fluctuation-dissipation relation). The pulse sequence 
with time interval ^ r cancels out the noise at frequen- 
cies lower than 1/r (dynamical decoupling). Therefore, 
the noise spectrum around the frequency 1/r dominantly 
contributes to the coherence time in the presence of the 
pulse sequence. 

Note that the infinite t limit of Eq. (fT3|) is consistent 
with Eq. dH): 



1/1 00 - 

hm = —^(0) y — 

^ 1=0 ^ 



ls{0) = ^, 



(15) 



where we use liniT-^oo ^^2;+! = 0, VZ in the first line, and 
X]/^o 1/(2^ + 1)^ = TT^/S in the second line. This equation 
enables us to interpret T2 as a generalization of Tj^^ 
into non-zero frequencies, which is the reason we call 
"generalized coherence time." 
Similarly, the infinitesimal t limit of Eq. (1131) is given 

by 



lim 



+0 



lim S{uj). 



(16) 



Measurement of noise spectrum. — Using the approx- 
imate relation (jl4p we evaluate the noise spectrum as 
follows. For a fixed value of 2t (inter-pulse time), we 
evaluate T2 by applying a sufficiently large number of 
pulses and measuring the asymptotic behavior of the co- 
herence. We repeat this procedure for other fixed values 
of 2t. Then, by plotting I/T2 against 7r/2r, we obtain 
5 as a function of w. 

Better evaluation of the spectrum is possible by us- 
ing the precise relation (jl3|) . For example, a functional 
form S^^'{uj) of the noise spectrum (which includes some 
fitting parameters) based on the above approximate eval- 
uation is phenomenologically introduced. Then we create 
a function F{tt/2t) similar to that on the right-hand side 
of Eq. (jl3p by summing the phenomenological function 
S^*" up to an appropriate cutoff L {L > 10 would be 
sufficient): 



F(7r/2T) 



±y. 



1 



(2^ + 1)^ 



5*^*((2/ + l)7r/2r). (17) 
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We finally fit F{tt/2t) to the experimental results {I/T2 
vs tt/2t) to obtain the values of the parameters for S^*". 

Finally, we estimate the frequency range for this 
method of noise measurement. To measure experimen- 
tally the asymptotic behavior of the coherence, we should 
apply a large number of pulses before the amplitude be- 
comes too small to detect. Hence, the inter-pulse time 
2t must be smaller than the 1/e decay time T2 of the 
coherence (in the presence of the SE pulse sequence). 
This implies that the lower bound on the frequency 7r/2r 
should be n/T2. The upper bound is determined by the 
shortest inter-pulse time that is experimentally available. 
The value of this time can be of the same order as that 
of the pulse duration time Tp. Hence, the upper bound 
on the frequency 7r/2r should be tt/t^. 

Concluding remarks. — In summary, we have described 
a method for obtaining the dephasing noise spectrum. 
The method is simple in the sense that we have only 
to apply sequences of equidistant 7r-pulses to the qubit 
(spin). The generalized coherence time, evaluated from 
the asymptotic exponential decay of the coherence in the 
presence of a sufficiently large number of pulses, is used 
for evaluating the spectrum. This method is applicable 
to a system interacting with several independent noise 
sources. In this case, the total spectrum is the sum of 
the individual noise spectra. 

This method is expected to be valid in a wide range 
of systems because we have derived it for three mod- 
els. The single-qubit noise spectrum plays an impor- 
tant role even in multi-qubit systems because it signif- 
icantly contributes to the dynamics of the systems (this 
is clearly seen when analyzing with the quantum master 
equation) . Extension of the present study remains a the- 
oretical task. We should analyze a model in which the 
system-environment coupling is given in a general form 
and/or in which there are both energy and phase relax- 
ations. The projection operator formalism 31[ would be 
helpful in analyzing a general model. 

This method should provide new insights into NMR 
experiments in condensed matter physics. So far the lon- 
gitudinal relaxation time Ti has been successfully used 
for investigating properties of electrons and nuclear spins 
in condensed matter. The method presented here should 
enable the use of the generalized coherence time T2 for 
investigating (other) properties of them. This is because 
the noise spectrum (evaluated from T2) includes infor- 
mation on the environment, the physical origin of which 
is electrons, other nuclear spins, and so on. In order 
to make this application useful, it is important to cap- 
ture the characteristic structure of the spectrum by ana- 
lyzing microscopic models that include interactions with 
nuclear spins (e.g., Fermi contact hyperfine and dipo- 
lar couplings). An experimental demonstration of the 
present method in NMR will be reported elsewhere 32 1. 
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SUPPLEMENTAL MATERIAL 

In this supplement, we give detailed derivations of 
Eq. (jlj for the three models for completeness. The essen- 
tials of the derivations are the same as those in Refs. [§- 

13, m. 



Classical-noise model 



The Hamiltonian of the classical noise model is 



(A.l) 



Here, 9k = n{tk+i-tk) + J*^'^' dt'^{t'). At time tk+i, an 
ideal 7r-pulse is applied about the x (y) axis, so 

Ps{tk+1 + 0) = D^(^y){±TT)ps{tk+l - 0)D.j;ly){±TT)\ 

(A.4) 

where the ±7r-rotation operator about the x (y) axis is 
Substituting Eq. (|X3l) into Eq. (|X4| . we obtain 

{+\'ps{tk+i + 0)1-) = ±{-\Utk + 0)\+)e^'\ (A.5) 
HUtk+i + 0)1+) = ±{+\kitk + 0)|-)e-^«'=. (A.6) 



where ^(t) is a stationary stochastic process with zero 

mean. In this model the noise spectrum defined by x (y) axis. By solving these recurrence relations, we 

Eq. (2)) can be rewritten as 



dt 



'mm) 



bath ■ 



(A.2) 



Here, the upper (lower) signs are for the rotation about 
the X (; 
obtain 

{ + \Mt^ + o)\-) 



The relation between the states at times tfc + (im- 
mediately after the fcth pulse) and tk+i — (immediately 
before the (fc + l)th pulse) is 



= (-l)"''(±|;5,(0)|T)exp 



n-l 



k=0 



where the upper (lower) signs are for even (odd) n, and Uy 
is the total number of rotations about the y axis. Then, 
Ps{tk+i — 0) = U{tk+i,tk)ps{tk + 0)tj{tk+i,tk)'' , (A. 3) from the state at time t (at which the signal is measured) 



where the unitary time evolution operator U is given by 



U{tk+i,tk) = exp 



At'H{t') 



Ps{t) = U{t,tn)ps{tn + 0)U{t,tny, we obtain 

{ + \Mt)\-) 



exp 1 -^(JzOk 



= (-l)"»(±|^.(0)|T)exp 



k=0 



Averaging over the stochastic variables, we obtain 



(-l)"V.,±T(0)exp 



ir!^(-i)^-+"(tfc+i -ife) 



fc=0 



exp 



k=o 



The normalized coherence W{t) is then given by (note that \psA I = \Ps. hi) 



W{t) = 



exp 



exp 



n 



k=0 
n+1 



rtk+i 






^bath 







dt'mftd') 



bath 



where ft{t') is defined by Eq. ([6]). 

If ^(i) is a stationary Gaussian process we obtain the desired result: 



W{t) 



exp 



exp 



di'i / dt^(C(t'i)e(t^))bath/t(ti)/t(i^) 

Aijj 
2^ 



bath 



(A.7) 



(A., 
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Analysis of pure dephasing model with quantum transformation, we obtain effective qubit operators and 

bath effective time evolution operators: 

Here we analyze a pure dephasing model described by 
the Hamiltonian ([T]). We assume that the Hamiltonian 

can be decoupled into the qubit and bath parts by a ^^(y) — ^^'^ ^x{y)^^"^'^ 



unitary transformation: 



e 



= e2*^^a,(,) = a,(,)e-2*^^, (A.IO) 

= ^<7,+7?^,th + (const.), (A.9) =exp[-iH''{tk+i-tk)/h]. (A.ll) 

where K and -ff^j^th ^'^^ respectively anti-Hermitian and 
Hermitian operators of the bath. Using the same unitary 

I 

The relation between the states at times tk + (immediately after the fcth pulse) and tk+i — (immediately before 
the (fc + l)th pulse) is 

p(tk+i-0)^eiip[-iH{tk+i-tk)/h]p{tk+0)ejcp[iH{tk+i-tk)/h]. (A.12) 
At time t/c+i, an ideal 7r-pulse is applied about the x (y) axis, so 

Kh+i + 0) = ^,(j,)(±7r)p(tfe+i - 0)£),(^)(±^)t, (A.13) 

where the ±7r-rotation operator about the x (y) axis is D^(^y-^{±T:) = ^ia^(^yy 
Substituting Eq. (|A.12|) into Eq. (|A.13|) . we obtain 

p{tk+i + 0) = (T^(y) cxp [-iH{tk+i - tk)/fi\ p{tk + 0) exp [iH{tk+i - tk)/h] ^^(y) , 

p'{tk+i + 0) = + 0)e-*^^ 

-e2*^*a,(,)C/'=(tfc+i,tfc)p^(tfe+0)?7^(ife+i,ifc)ta,(^)e-2-^*, 

{+\p{tk+i+0)\-) = (+|e-*^*p^(tfc+i +0)e*^^|-) 

^ ±e'^^('^+^~'-)f^-f +{^\p{tk + 0)|+)ffcY^+t, (A.14) 
{-\p{tk+i+0)\+) = ±e~'^<-''<'+'-*''^e-^ e-'^^-''^-^^'^*''+'-*'''^/''e^ 

= ±e-^(*'=+-*'=)f^;\':,-(+|p(t, +0)|-)f^f/-t. (A.15) 

Here, f^f^+ = e-'Hi^^i,itk+i-tk)/h^-k ^ Tfc+iV = e--^e-*-^bath(*<=+i-*<e)/'ie^; the upper (lower) signs are for 
rotation about the x (y) axis. 

By solving these recurrence relations (jA.14p and (jA.151) . we obtain 

(+|p(t„ + 0)1-) =(-l)nyp-^I:'^=i(-l)'+"n{t, + l't,)fh.t^^ . . . j,bathT 

X (±|p(0)|T)rl^r^t . . . rbath-_^tf ba^th+t^ (A.16) 

where the upper (lower) signs are for even (odd) n, and Uy is the total number of rotations about the y axis. From 
the state at time t (when the signal is measured) p{t) = exp[—iH{t — tn)/h\p{tn + 0) exp [iff (t — tn)/h\ , we have 

(+I/5WI-) = (+1 cxp[-^ff(^ - t„)/h] \+){+\pitn + 0)|-)(-| cxp[-^ff(^ - t„)/h] I-) 

= e-»^^(*-*")e-i<'e-^^^L«,(*-t.)/'^e*(+|^(t„ + 0)|-)e^"e*^-'^(*-*")/''e-^, (A.17) 

where we have used {-\exp[-iH{t-tn)/h]\+) = {+\cxp[-iH{t-tn)/h]\~) =0. Combining Eqs. (|A.16p and (|A.17|) 
and taking the partial trace of the bath, we obtain 

p+-{t) Trbath(+|/o(t)|-) 

= (-l)"«e~'^^=o(-i)'*""(*'=+i-*'')Trbath f(±|/5(0)|T)e("i)"2^1 , (A.18) 
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where to = 0, tn+i = t, and 



L = k + ^(-i)^-2X(tfc) + {-iY+^k{t). 



(A.19) 



fe=0 



Here, k{T) = e'"^.a,T/hj;^^-tH^^^^^T/h^ jf ^^le initial state is a decoupled one, p(0) = Ps(0)Pblth' 



^oq (-1)"2L 



P±t(0), 



(A.20) 
(A.21) 



Spin-boson model 



Spin-spin bath model 



The Hamiltonian of the spin-boson model is 



i/bath, (A.22) 



where = J2j + ^j) ^'^'^ ^^bath = J2j ^'^j^j^j- The 
spectral density of the bath J{ijo) is defined by 



(A.23) 



for oj > 0. The relation between J(w) and S'(w) is 

S{Lo) = TiJ{\oj\){2rn,{\u:\) + l). (A.24) 

Here, nb(a;) = l/(e^'''^ — 1) is the boson occupation fac- 
tor. 

The spin-boson model satisfies condition (jA.9|) if we 

= ^a^+HU, + AE, (A.25) 

where ff^ath = -f^bath Z]j^j^]^jj and Ai? = 
— dwJ(aj)/4w. In this case, Eq. (|A.19p becomes 

A 



The Hamiltonian of the spin-spin bath model is 

h 



H 



k 



bath 1 



(A.27) 



where | = J2j + ■s'^) and 7?bath = Y.j{^/'^)^3^1- 

The spectral density of the bath J{u)) is defined by 

JH=^M-'^(^-^,), (A.28) 
i 

for w > 0. The relation between J(w) and is 

S'(w) = 47rJ(|a;|). (A.29) 



The spin-spin bath model satisfies condition (IA.9I) if 
we set k = (z/2) Y,j "jS^: 



-cf, + H, 



bath' 



(A.30) 



where H^^^^^ = (ft/2) ^^.(wj cosa^ + sin , and 
aj = arctan(2^j7a;j). In this case, Eq. (|A.19I) becomes 



siReA(7,)-s^,Im/,(7,) , (A.31) 



^ = ^ [/t('^j)^] - /t*K>j > (A.26) where 7^ = yw|T4^, and /t(w) is given by Eq. ©. 



where /t(w) is given by Eq. ([5|). Then, by applying the 
Wick theorem, we transform Eq. (jA.21[) into 



' bath 



W{t) = exp(2(L2) 

dc^JH|/iH|'(2nbH + l 



exp 



Using Eq. (jA.24l) , we finally obtain the desired result (jH) . 



Then 

j 

where 6^ = aj7j|/t(7j)|, = Re/t(7j)/|/t(7j)|, and 
Y,- = -Im/t(7j)/|/f(7j)|. Therefore Eq. (|A.2ip becomes 

W{t) = J||cos6ij| 

3 

= exp ^ In I cos 6j\ \ . (A.32) 



If fij <C ujj (the coupling between the qubit and the 
bath is weak), aj = 2fij/ujj + 0(^{iij /cjj)'^) , 7j = uij + 
0{{fi,/u,f), and /t(7,) = ft{oj,) + 0{{fi,/uj,)^). There- 
fore, 9j = 2fj,j\ft{ujj)\ + 0{^{^j /ujj)^) , and ln|cos0j| — 
— 2^||/t(cjj)|^ + 0(^{^j/ujjY). Substituting these equa- 
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tions into Eq. (jA.32l) . we obtain 

= exp(^-2^ dujJ{uj)\~ft{uj)\^ + 0{{^iJ|uJ,f)^ . 
Using Eq. (|A.29p , we finally obtain the desired result (jH) . 



